Exact solution of the Percus— Yevick integral equation for fluid mixtures of hard 

hyperspheres 
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Structural and thermodynamic properties of multicomponent hard-sphere fluids at odd dimen- 
sions have recently been derived in the framework of the rational function approximation (RFA) 
[Rohrmann and Santos, Phys. Rev. E 83, 011201 (2011)]. It is demonstrated here that the RFA 
technique yields the exact solution of the Percus- Yevick (PY) closure to the Ornstein-Zernike (OZ) 
equation for binary mixtures at arbitrary odd dimensions. The proof relies mainly on the Fourier 
transforms Cij(fc) of the direct correlation functions defined by the OZ relation. From the analysis 
of the poles of Cij(fc) we show that the direct correlation functions evaluated by the RFA method 
vanish outside the hard core, as required by the PY theory. 
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I. INTRODUCTION 

The understanding and interpretation of thermody- 
namic properties and the structure of real dense gases 
and liquids rely mostly on the use of models of fluids and 
the application of approximate theories [ij. One of the 
most successful and widely used approximate methods is 
the Percus- Yevick (PY) theory. It is an integral equa- 
tion theory based on the Ornstein-Zernike (OZ) equa- 
tion m coupled with the PY closure Q . Systems of rigid 
nonattracting hard spheres (HS), as well as their modified 
forms (e.g., square- well and sticky particles), represent 
useful fluid models for the study of liquids and colloidal 
systems In this context, it is noteworthy that exact 
solutions of the PY integral equation were early obtained 
for pure 0, Q and multicomponent HS fluids. 

Recently, we have shown that the thermodynamic 
and structural properties of single and (additive) multi- 
component fluids of hard hyperspheres at odd dimensions 
can be studied by means of an analytical technique, the 
so-called rational function approximation (RFA). This 
method is formulated in terms of the Laplace transform 
of a polynomial times the radial distribution function and 
leads directly to a system of algebraic equations that have 
analytical solution at low spatial dimensionality (odd 
d < 7 for pure fluids and odd c? < 3 for mixtures) and 
must be solved numerically at higher odd dimensions. 

In the single-component case, it was shown that, in 
its simplest formulation, the RFA method recovers the 
known solution of the PY closure to the OZ equation Q. 
In the multicomponent case, the equivalence between the 
solutions of the two approximations has been shown for 



binary mixtures in three dimensions [To| and five dimen- 
sions 0] . Although compelling arguments were presented 
in Ref. Q in favor of the extension of the equivalence 
to general odd d, a rigorous proof was not attempted. 
Therefore, whether or not the PY integral equation is ex- 
actly solvable for mixtures of additive hard hyperspheres 
with arbitrary d ~ odd (as well as the structure of the 
solution in the affirmative case) remained, to the best of 
our knowledge, an open problem. 

The purpose of the present paper is to give a formal 
proof that the PY and RFA solutions are indeed equiv- 
alent for binary mixtures in arbitrary odd dimensions. 
The proof is based on the analysis of the Fourier trans- 
form Cij (fc) of the direct correlation function (r) . which 
is defined by the OZ equation and is a key quantity to 
be determined by both approximations. The equivalence 
between RFA and PY theories is ensured by first deriving 
the functional structure of Cij{k) as given by the RFA. 
The direct correlation function Cij(r) is then calculated 
by the theory of residues, showing that it is strictly short 
ranged, as expressed by the PY closure. 

The organization of the paper is as follows. Section 
im provides a description of the fiuid structure equations 
and basic definitions of the two approximations. Section 
mil particularizes the RFA scheme to binary mixtures of 
HS at odd dimensionalities. Section IIVI deals with the 
analysis of the direct correlation function and its short- 
range behavior as derived from the RFA theory. We close 
the paper in Sec. |V] with some concluding remarks. 



II. FORMULATION 



A. Basic quantities 
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Let US consider an A/'-component fluid mixture of ad- 
ditive c?-dimensional HS. Let {xi;i = 1,...,A^} be the 
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set of mole fractions, {ai;i = 1,...,A^} be the set of 
diameters, and p be the total number density. 

The radial distribution function gij(r), the total cor- 
relation function hij{r), and the direct correlation func- 
tion Cij{r) corresponding to particle pairs of species i and 
j are the primary linkage between thermodynamic and 
structural properties and particle interactions of multi- 
component fluid mixtures. They are not independent 
quantities. In particular hijir) = gij{r) — 1, while Cy(r) 
is defined in terms of hij(r) by the OZ equation, which is 
conveniently written in Fourier space and in matrix form 
as 



(ii) Short wave-number expansion: 



G, 



_ (d-2)!! V" rr ™ 



where 



m=0 



(2.7) 



(2.8) 



Zik) 



I 



U{k)- l + h(fc) 



(2.1) 



Here, I is the J\f x J\f unit matrix and c{k) and h{k) 
are Af x Af matrices with elements p^XiXjCij{k) and 

Py^XiXjhij{k), respectively, where the hat symbol (^) de- 
notes the Fourier transform. 

The Fourier transform of the total correlation function 
can be evaluated from 



h,j{k) 



i-2n 



s—ik 



,(d-l)/2 



(2.2) 



where i is the imaginary unit and Gij{s) is a Laplace 
functional of the radial distribution function defined by 



Gij{s) = / drrgij{r)0n{sr)e 
Jo 



with 



(2.3) 



(2.4) 



dr hij(r)r'" 



(iii) Low-density expansion: 



G.j{s) 



=1 

2^ 



(2.9) 



+ p'^xi / Arr9n{sr) 



(2.10) 



where VLa^bif) is the intersection volume of two hy- 
perspheres of radii a and h whose centers are sepa- 
rated by a distance r < a + 6. 

Conditions (i)-(iii) are consequences of the hard-core 
interaction and, therefore, are directly related to the 
step-function structure of gij[r) [Eq. p.Sp ]. Because 
oin,2q+i = for (7 = 0, ...,n— 1 3, condition (ii) guaran- 
tees through (|2.2p that the correlation functions remain 
bounded at the limit of zero wave-number. 



%(0) < oo. 



(2.11) 



being the reverse Bessel polynomial of degree n = [d — 
3)/2. The function Gij{s) is a fundamental quantity in 
the sense that it contains all the information about the 
thermodynamic properties of the fluid and its knowledge 
gives structural information equivalent to pair radial dis- 
tributions, pair correlation functions, and structure fac- 
tors. 

Since the HS particles are not allowed to penetrate each 
other, the radial distribution function vanishes inside the 



core. I.e. 



(2.5) 



where Uy- = ^{ai + aj ) is the contact distance for hyper- 
spheres of species i and j. In addition, the exact Laplace 
functionals Gij{s) for d-dimensional HS satisfy the fol- 
lowing asymptotic relations 0]. 



(i) Long wave-number limit: 



lim s 

s— >oo 



(5-d)/2^ 



»G,,(s) 



_ Jd-l)/2 



5u(4)- (2-6) 



and this in turn assures that the isothermal compress- 
ibility of the fluid takes finite values @. 



B. The PY approximation 

The PY approximation to a classical fluid is obtained 
from the OZ equation (|2.ip supplemented by a particu- 
lar closure relation between gij{r) and Cij{r). For hard- 
hypersphere systems, the PY closure reduces to Eq. ()2.5|) 
and the assumption that the direct correlation function 
is short ranged, namely, 



Cijir) =0, r > ai 



(2.12) 



Therefore, if an approximate radial distribution function 
gij (r) satisfies the hard-core condition (|2.5p and its asso- 
ciated direct correlation function, as obtained from Eq. 
(|2.ip . verifies the condition (|2.12[) . it is necessarily a so- 
lution to the PY theory. 
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C. The RFA 

In contrast to the PY approach, the RFA method is 
based on the Laplace functionals Gij{s) rather than on 
the OZ equation. Specificahy, the RFA provides an an- 
alytical representation of Gij{s) that complies with the 
consistency conditions (i)-(iii) Q. This analytical ap- 
proximation reads 



(2.13) 



where L(s) and B(s) are J\f x Af matrices with elements 



n+l 



m=0 



s di 



VpiE Lij{s) + Pij{s). 



(2.14) 



(2.15) 



In Eq. (|2.15p , (5^ is the Kronecker delta symbol, pi = pxi 
is the partial number density of species i, and 



n+l 
?n— 



d—m 

E 

.i=0 



i-<y^sY 
ei 



(2.16) 



Note that the matrices L(s) and P(s) have polynomial 
dependencies on s of degrees n+l = (d— 1)/2 and 2n+3 = 
d, respectively 

The coefficients L^™\ m = 0, l,...,n + 1 = (d - l)/2, 
may depend on the fluid density, the particle diameters, 
and the component abundances, but they are indepen- 
dent of s. Those coefficients are determined from a set 
of 71 -|- 2 algebraic matrix equations stemming from the 
requirement of condition (|2.7p . The reader is urged to 
consult Ref. Q for further details. 

Note that Eq. (|2.13p can be rewritten as 



G.J is) 



F.J is) 



where 



F(s)EE L(s).adj[B(s)], ^(s) = |B(s) 



(2.17) 



(2.18) 



Here, adj(A) and |A| refer to the adjoint and the deter- 
minant, respectively, of a matrix A. Once the coefficients 
Lj-™^ are obtained as functions of p, {ci}, and {xi}, the 
functionals G,;j(s) are fully determined from Eqs. (|2.13p 
or p.l7p . Then, the total correlation functions in k space 
are obtained from Eq. (|2.2p , which can be rewritten using 
(PTT]) as 



Kj{k) 



Ri] (ifc) 

' m) Pi-iky 



(2.19) 



where 



i?y(s) = F,j(s)/3(-s)e-"-^ + F„(-s)/3(s)e"-^ (2.20) 



Finally, to obtain the direct correlation functions in the 
configuration space, one makes use of Eq. (|2.ip and the 
inverse Fourier transform 



(27r 



|-(d+l)/2 poo 



„d-2 



dk kcij ik)6n iikr)e' 



-ikr 



(2.21) 

It is worth noting that, as shown in Ref. (9[, the Laplace 
functionals obtained from (|2.13p satisfy the asymptotic 
behaviors given by Eqs. p.6p - p.l0p . In particular, the 
physical hard-core requirement given by Eq. (j2.5p is veri- 
fied. Therefore, in order to prove the equivalence between 
the PY and RFA approaches it will be sufficient to show 
that the direct correlation functions calculated with Eqs. 
(|2.ip . (|2.19p . and (|2.2ip are short ranged, as required by 

Eq. dug). 

Before closing this section, let us derive some proper- 
ties for small s that will be useful in Sec. IIVI By ex- 
panding the exponential in Eq. (j2.15p and inserting Eq. 
one gets 



n+l 



Therefore, i?y (s) 
nant is 



vp, 

m=0 

s^lS,. 



i~<J.s)' 



T (™) „" 



E 

'=d-m+l 

(2.22) 

0(5)], SO that the determi- 



/3(s)-s^'^[l + 0(s)]. 
Next, from Eqs. dUT]) and (PTT)) we get 



(2.23) 



e-'^--^^., (.) = id- 2)!!^ [1 + s^F^.is^) + Ois'')] , 

* (2.24) 
where J^ijis^) is a polynomial of degree n = (d — 3)/2 in 
whose explicit form will not be needed here. In fact, 
the polynomial Tijis'^) cancels in Eq. p.20p . resulting 
in %(s) = ^(s^Vd 

obtains 



From Eqs. (^^^ and (p:^ one 



F,,(.s) = (d-2)!!s(^-i)'^[l + 0(.s)] 



(2.25) 



Particularizing to the binary case (A/" = 2) and taking 
the determinant on both sides of Eq. (|2.24p , we get 



-2cri2S I 



Hs)\ 



[j-(s") + Ois 



d-2\ 



(2.26) 



where J-^(s^) is a polynomial of degree n = (d — 3)/2 in 
s^. Combination of Eqs. ([^^ and (P?^ yields 

e(2^.-^.>|L(-.)|F.,(.) = ^^^^ [£,^.(,2) 

+Ois'^-^)] , (2.27) 

where Cijis^) is again a polynomial of degree n = (d — 
3)/2 and use has been made of the property 

|F(s)|=/3(s)|L(s)|, (2.28) 
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which follows from Eq. (|2.18|) . Finally, Eq. (|2.27p implies 



(2.29) 



where in the last equality we have taken into account Eq. 

Km . 



III. BINARY MIXTURES 

Henceforth, we consider an additive system of two com- 
ponents {Af = 2). From Eq. (12. ip . we have 

^ feii(fc)+P2|h(fc)| ^ hi2{k) 
cii(k) = ^^77;^ , ci2(fc) = , (3.1) 



Dik) 



D{k) 



with 



D{k) = 1 + pi/iii(fc) + p2h22{k) + piP2Hk)\- (3.2) 

The expressions for C22(fc) and C2i(fc) are obtained from 
those of cii(fc) and Ci2(fc) by exchanging the subscripts 
1 and 2. 

Equations (|3.ip and (|3.2p are general and valid for any 
binary mixture. Now we particularize to the RFA. Using 
Eq. (|2.19p . one obtains 



D{k) 



V{ik) 



/3(ifc)/3(-ifc) ' 



(3.3) 



c,,(/0 ^ ,. -^"(^^') + '^P2|R(ifc)|/[/3(ifc)/3(~ifc)] ^ ^3_4^ 



2?(iA;) 



Ci2(fc) = 1/ 



i?i2(ifc) 

I?(ifc) ' 



(3.5) 



where 



2?(s) = /?(s)/3(-s) + z^[pii?ii(s) + p2i?22(s)] 



|R(s) 

-I^ /91/92- 



■/?(s)/3(-.)- 
Explicit calculation using Eq. (|2.20p gives 

|R(s)l 



$(s) + $(-s) 



(3.6) 



(3.7) 



V{s) = P{s)l3i-s) + A{s) + A{~s), (3.5 



where 



$(5) = |L(s)|/3(-s)e-(-i+-^)^ - Fi2(s)F2i(-s) 

+Fn(s)F22(-s)e-(-l-^-)^ (3.9) 



A(5) = z^/3(-5)[piFn(s)e-'^i^ +p2F22(5)e-'^=1 

WpiP2Hs). (3.10) 



In Eq. dm we have used Eq. (H^S]). 

Thus far, we have not used in this section the explicit 
form of the matrix B(s), Eq. (|2.15p . Evaluation of Eq. 
(|2.18p for binary mixtures yields 



where 



Fnis) - iyp2\l{s)\e-''-' ^ Fi2(s) 

F2i{s) F22is) - i^pi\Lis)\e--^-' 

(3.11) 



Fn(s) = Ln(s)[s'' + P22(s)] - ^12(5)^21(5), (3.12) 

Fi2(s) = Li2{s)[s'' + Puis)] - iii(s)Pi2(s). (3.13) 

Of course, ^22(5) and -^21(3) are obtained from Eqs. 
p.l2p and p.l3p by exchanging the subscripts 1 and 2. 
In addition, the determinant /3(s) of the matrix B(s) can 
be expressed as 

/3(s) =^(s)+z.Vp2|L(s)|e-(^i+'^^)^ 

-iypihi{s)e'^'' - iyp2F22{s)e-^'', (3.14) 



with 



Pis) EE S'^ + S^[Pu{s) + P22(.S)] + |P(.S)|. 



(3.15) 



We can observe that -Fij(s) and /3(s) are polynomials 
in s of degrees 3n + 4 = (3d - 1) /2 and 2(2n + 3) = 2d, 
respectively. In contrast, the functions i^ii(s), ^22(5), 
Pis), $(s), and A(s) contain exponential terms. The 
key point, however, is that those exponential terms com- 
pensate exactly in the function I?(s). Using Eqs. p.9p . 
(|XTU1) . ([XTT|) . and ([XTi)) in Eq. (g^ one finds, after 
some algebra. 



Vis)^l3is)Pi~s) + :.^pipi\Lis)\\Li-s)\ 
2 

-1/2 p,PjF,jis)Fj,i-s). 



(3.16) 



Thus, X'(s) is an even polynomial of degree 4(2ri -I- 3) = 
4d. 

The direct correlation functions in Fourier space are 
given by Eq s. ([ Mt and (|331) . From Eqs. (E^O]), ([SJ]) . 
(1331), ((XTT|) . and (|XTi)) it is possible to get 



Cij(fc) 



X'(ifc) 



7',j(ifc)e'"-'^+7'y(-iA:)e-'"-'= 



~Q^,i^k)e 



i{(Ti —(Tj)k/2 



Q.,i~\k)^ 



{(jj 



,)fc/2 

(3.17) 



where 'Py(s) and Sy(s) are polynomials of degrees 7n + 
10 = (7d - l)/2 and 2(3n -f 4) = 3d - 1, respectively, 
given by 

^ii(s) = Fui~s)Pis) - y''pl\Li-s)\F22is), (3.18) 
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Vi2{s) = Fi2(-s)/3(.s) + i/Vp2|L(-s)|Fi2(s), (3.19) 

Qii(s) = -V piFii{-s)Fii{s) + P2Fi2{-~s)F2i{s) 

WpwlUs)\U~s)l (3.20) 

Qi2{s) = -v piFn{-s)Fr2{s)+P2Fi2{-s)F22{s) ■ 

(3.21) 

Again, 7^22(5), 7^21 (s), Q22(s), and Q2i(s) are obtained 
by exchanging the subscripts 1 and 2. 

IV. ANALYSIS 

In the derivation of the results in Sec. IIIII we have not 
needed to use cither Eq. (|2.16p or the conditions that the 

coefficients l[J^^ must satisfy [which are summarized by 
Eq. ([2?24| ]. In fact, Eq. ([STT]) alone is not sufficient to 
prove the PY condition ([2TT2|) . 

As said before, the RFA method guarantees that hij{k) 
takes finite values at the limit of zero wave-number, Eq. 
(PTT|) . Thus, from Eq. D{k) remains bounded at 

k — Q and, according to Eq. (|3.3p . we have 



1- ^(^) 
™o /3(s)/3(-s) 



= £1(0) < 00. 



(4.1) 



Taking into account that, in the binary case [M = 2), 
j3{s) = s^^ for small s [sec Eq. ^(272^ ]. we have V{s) = 
C(s^''). Since, according to Eq. p.l6p . I?(s) is a poly- 
nomial of degree 4(i, we conclude that the I?(s) is just a 
pure power law, i.e.. 



V{s) = D(0) 



4d 



(4.2) 



This is the crucial result allowing one to prove Eq. (j2.12p 
from Eq. p.l7p . Before proceeding to the proof, let us 
first simplify Eq. (|3.17p a little more. 

As shown by Eqs. ((XT5|) " ((??^ . Vij{s) and Q,j(s) arc 
polynomials of degrees 7n -f 10 = {Id — l)/2 and 2(3n + 
4) = 3(i — 1, respectively. Going back to the quantities 
Fij{s) and /3(s), Eqs. (|37l8l) - (|33T|) can be rewritten as 

2 

"Prjis) = F,j{-s)P{s) + vF,j{-s)Y,Pke-'"''Fkk{s) 



k=l 



+v piP2e' 



+e-(2-i=-'^-)^|L(s)|F„(-s)' 



(4.3) 



-J^5]pfcF»fc(-.s)Ffc,(,s) 
fe=i 

-z.ViP2e-("'^"^'-^/' [e(2----->|L(-s)|f^,,(s) 
+e-(2-i=-"-'^|L(s)|F,,(-s)l . (4.4) 



While Eqs. (|4.3p and (j4.4p conceal the polynomial char- 
acter of 'Pij(s) and Qij{s), they show, with the help of 
Eqs. ([2:231 . ([225]), and (p^ . that 



7'.y(,s) = 0(s^''), Q,,(s) = 0(s^'') 



(4.5) 



This result, along with Eq. (|4.2p . allows us to rewrite Eq. 
([STf]) as 



%(^^) = i [7'.,(iA:)e--'= -H7',,(-ifc)e 



+ Q,j(ifc)e'("--"^)''^/2 + Q,j(-iA:)e'('^^-"')*'-/2] , 

(4.6) 



where 



^y(5) 



D{0) 



is a polynomial of degree 3n + 4 = 3(o?— l)/2 and 



(4.7) 



(4.8) 



is a polynomial of degree 2(n-|-l) = d— 1. 

Equation (|4.6p is the main result of this paper. It 
provides the functional structure of the Fourier trans- 
form Cij{k) in the RFA approach. The direct correla- 
tions function in the configuration space Cij{r) can be 
calculated by application of the residue theorem com- 
bining Eqs. ^(TH) and (jM]). Since %(0) = finite, the 
integrand in Eq. p.2ip is regular along the real axis and 
so we can distort the integration path in the complex k 
plane by going around the point A; = from below. Next, 
the integral in Eq. (|2.2ip decomposes into four contribu- 
tions with integrands headed by e~^^'^'^~'^^'\ Q-^Kr+iyij) ^ 
g-ifc[r+(^,-a,)/2]^ and e-"=['-+('^'-'^^)/2l, respectively, with 
each integrand having a single pole at fc = of order 
2d—l. li < r < aij, the first integral must be closed 
with an upper half circle of infinite radius, and the residue 
theorem yields a nonzero value. If < r < [dj —(7^1/ 2, an 
additional nonzero contribution results from the third or 
fourth integral, depending on whether cr, > aj or aj > Oi , 
respectively. On the other hand, if r > aij , wc must close 
the path with a lower half circle. As a consequence, the 
four contributions vanish, and so the RFA method yields 

Eq. dang). 

This completes the proof on the equivalence between 
the PY and RFA solutions for binary mixtures of hard 
hypersphcrcs at odd dimensional space. 



V. CONCLUDING REMARKS 

PY and RFA theories are, in principle, alternative 
methods for calculating thermodynamic and structural 
functions of HS systems. The PY theory consists of the 
OZ relation (|2.ip supplemented with the hard-core con- 
dition (j2.5p and the genuine PY closure (|2.12p . In the 
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RFA approach, however, one proposes a specific form, 
Eqs. (|2.13p - (|2.16p . for the s-dependence of the Laplace 
functional Gij{s) defined by Eq. p.3|) . This specific form 

includes n + 2 = (d + l)/2 coefficients, which are 

determined, in consistency with Eq. (j2.7p . by requiring 
the independent term in the Taylor series expansion of 
s^Gij{s) to be {d - 2)!! and all the coefficients of s^-J+i 
with g = 0, . . . , n to vanish. 

In this paper we have shown by a direct verification 
that both methods arc fully equivalent for binary mix- 
tures of additive hard hyperspheres at odd dimensions. 
The proof is based on the analysis of the Fourier trans- 
form Cij{k) of the direct correlation function Cy (r) and 
proceeds along two main stages. In the first stage, use 
of Eqs. (PJ^ ~ (PT^ has allowed us to derive Eq. pTTTl) . 
where Vij{s), Qy(s), and 2?(s) are polynomials. This 
result applies regardless of the values of the n + 2 coef- 
ficients In the second stage, we have proved that 
enforcement of Eq. ([2J)) implies Eqs. (|4?2|) . (|4J)) . and 
(g^, so that Eq. ((XT7| simplifies further to Eq. (gH). 
Application of the residue theorem then yields Eq. (|2.12p , 
which completes the proof. An interesting feature of the 
proof is that the explicit expressions for the coefficients 
L,-™'' as functions of the physical parameters of the fluid 
(density, mole fractions, and particle diameters) arc not 



needed. 

Although the proof presented in this paper has been 
restricted to the binary case (A/" = 2), we conjecture that 
the structure of Eq. (|4.6p , and hence the validity of Eq. 
(|2.12p . remains valid for any number Af of components 
in the framework of the RFA. 

It is worth mentioning that, as done for three- 
dimensional mixtures [loj and for d-dimensional one- 
component systems [1] , the RFA scheme can be extended 
beyond the FY level by adding extra terms i^""''^'' in Eqs. 
(|2.14[) and (|2.16p . and replacing Sij by (1 +us)Sij in Eq. 

(|2.15p . The Af^ + 1 free parameters and u can be 

fixed by imposing given expressions for the contact values 
gij{cr^j) and the thermodynamically consistent isother- 
mal compressibility. 
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